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Majorana fermions hold promise for quantum computation, because their non-Abelian braiding 
statistics allows for topologically protected operations on quantum information. Topological qubits 
can be constructed from pairs of well-separated Majoranas in networks of nanowires. The coupling 
to a superconducting charge qubit in a transmission line resonator (transmon) permits braiding of 
Majoranas by external variation of magnetic fluxes. We show that readout operations can also be 
fully flux-controlled, without requiring microscopic control over tunnel couplings. We identify the 
minimal circuit that can perform the initialization-braiding-measurement steps required to demon- 
strate non-Abelian statistics. We introduce the Random Access Majorana Memory, a scalable circuit 
that can perform a joint parity measurement on Majoranas belonging to a selection of topological 
qubits. Such multi-qubit measurements allow for the efficient creation of highly entangled states 
and simplify quantum error correction protocols by avoiding the need for ancilla qubits. 



After the first signatures [1-4] of Majorana bound 
states in superconducting nanowires [5-7], the quest 
for non-Abelian braiding statistics [8-11] has intensi- 
fied. Much interest towards Majorana fermions arises 
from their technological potential in fault-tolerant quan- 
tum computation [12-16]. Their non-Abelian exchange 
statistics would allow to perform quantum gates belong- 
ing to the Clifford group with extremely good accu- 
racy. Moreover, topological qubits encoded non-locally 
in well-separated Majorana bound states would be re- 
silient against many sources of decoherence. Even with- 
out the applications in quantum information processing, 
observing a new type of quantum statistics would be a 
milestone in the history of physics. 

The two central issues for the application of Majorana 
fermions are (i) how to unambiguously demonstrate their 
non-Abelian exchange statistics and (ii) how to exploit 
their full potential for quantum information processing. 
The first issue requires an elementary circuit that can 
perform three tasks: initialization of a qubit, braiding 
(exchange) of two Majoranas, and finally measurement 
(readout) of the qubit. In view of the second issue, this 
circuit should be scalable and serve as a first step towards 
universal fault-tolerant quantum computation. 

Here we present such a circuit, using a superconduct- 
ing charge qubit in a transmission line resonator {trans- 
mon [17-20]) to initialize, control, and measure the topo- 
logical qubit. In such a hybrid system, named top- 
transmon [21], the long-range Coulomb couplings of Ma- 
jorana fermions can be used to braid them and to read 
out their fermion parity [21, 22]. While there exist sev- 
eral proposals to braid or measure Majorana fermions in 
nanowires [11, 21-31], combining braiding and measure- 
ment without local adjustment of microscopic parame- 
ters remains a challenge. We show that full macroscopic 
control is possible if during the measurement one of the 
Majorana fermions is localized at a T-junction between 
three superconducting islands (see Fig. 1). All three 
steps of the braiding protocol, initialization-braiding- 
measurement, can then be performed by adjusting mag- 
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Figure 1. Two circuits that can demonstrate non-Abelian 
statistics, by the initialization, braiding, and measurement 
of pairs of Majorana bound states (circles). Braiding is per- 
formed twice to flip the fermion parity of 7^ and 7s [13]. 
Majoranas that can be coupled by Coulomb charging energy 
are connected by a thin line; the line is solid if the Majoranas 
are strongly coupled, and dashed if they are uncoupled. A 
thick line indicates tunnel coupling of Majoranas. The T- 
shaped circuit of Ref. [11] (left column) requires control over 
tunnel couplings, while the 7r-shaped circuit considered here 
(right column) does not, because both readout and braiding 
involve a Majorana localized at a T-junction. 



netic fluxes through split Josephson junctions. Because 
local control of microscopic parameters is not necessary, 
our scheme is less sensitive to problems arising from elec- 
trostatic disorder and screening of gate voltages by the 
superconductor. 

This design principle of flux-controlled braiding and 
measurements can be scaled up from a minimal braiding 
experiment setup to a multi-qubit register that supports 
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Figure 2. Panel a): Minimal circuit for flux-controlled demonstration of non-Abelian Majorana statistics. Two large supercon- 
ducting plates form a Cooper pair box in a transmission line resonator, i.e. a transmon qubit. Three smaller superconducting 
islands are embedded between the two transmon plates. Each superconducting island contains a nanowire supporting two 
Majorana bound states. At low energies, the three overlapping Majorana bound states at a T-junction form a single zero mode 
so that effectively the system hosts six Majorana bound states, labeled 7a, 7s, 7c, Id, "1e, and 7f. The Coulomb couplings 
between the Majorana fermions can be controlled with magnetic fluxes This hybrid device can measure the result of the 
braiding operation as a shift in the microwave resonance frequency when the fermion parity i'^A'ys switches between even and 
odd. Panel b): Sequence of variation of fluxes during the initialization (steps 0-2), braiding (steps 3-8) and measurement (step 
9). Panel c): Illustration of the steps required for initialization, braiding and measurement. To unambiguously demonstrate the 
non-Abelian nature of Majoranas, one needs to collect statistics of measurement outcomes when the adiabatic cycle describing 
the braiding operation (steps 3-8) is repeated n times between initialization and measurement. The probabilities of observing 
changes in the cavity's resonance frequency, Pflip, for different values of n should obey the predictions summarized in the table. 
The sequence of probabilities shown in the table repeats itself periodically for larger values of n. 



a universal set of quantum gates and allows measurement 
of any product of Pauli matrices belonging to a selection 
of topological qubits. Multi-qubit parity measurements 
are a powerful resource in quantum information process- 
ing, allowing for the efficient creation of long-range entan- 
glement and direct measurement of stabilizer operators 
(thus removing the overhead of ancilla qubits in quan- 
tum error correction schemes). Because the data stored 
in the register can be accessed in any random order, it 
truly represents a Random Access Majorana Memory. 

The structure of the paper is as follows. In Sec. I we 
present the circuit that can demonstrate the non-Abelian 
Majorana statistics. In Sec. II we take a longer-term 
perspective and describe the Random Access Majorana 
Memory, whose potential for quantum computation is 
discussed in Sec. III. Finally, we conclude in Sec. IV. 



For the benefit of the reader, we include more detailed 
derivations and discussions in the Appendices. 



I. MINIMAL CIRCUIT FOR THE 
DEMONSTRATION OF NON-ABELIAN 
STATISTICS 

To demonstrate non-Abelian Majorana statistics one 
needs to read out the parity of two Majoranas, and 
7b, and braid one of these Majoranas 7^ with another 
one, 7c- We seek a transmon circuit that can combine 
these operations in a fully flux-controlled way, by acting 
on the Coulomb coupling of the Majoranas. Since 75 
must be coupled first to one Majorana (for the braiding) 
and then to another (for the readout), it must be able 
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to contribute to two different charging energies. This is 
possible if 7b is locahzed at a T-junction between three 
superconducting islands. 

We thus arrive at the minimal circuit shown in Fig. 2a. 
It consists of five superconducting islands, each contain- 
ing a nanowire supporting two Majorana bound states, 
enclosed in a transmission line resonator. The two big- 
ger superconductors form a transmon qubit and the three 
smaller islands are embedded between the two transmon 
plates. The Josephson couplings between the islands can 
be controlled by magnetic fluxes {k = 0, 1, 2, 3). The 
nanowires form a 7r-shaped circuit, with two T-junctions 
where three Majorana bound states belonging to adjacent 
superconductors are tunnel-coupled. At low energies the 
three overlapping Majorana bound states at a T-junction 
form a single zero mode, so that effectively the system 
hosts six Majorana bound states, 7a, 7b, ■••,7f- 

The three relevant energy scales for the device are 
(i) the charging energy Ec,k = e^ /2Ck determined by 
the total capacitance Cu of the four upper superconduc- 
tors in Fig. 2a, (ii) the Josephson energies -Ej,fc($fc) = 
-£'j,/c(0) cos(e$fe//i), and (iii) the Majorana tunnel cou- 
plings i?M at both T-junctions. For strong Josephson 
coupling, E}^k 3> Ec^k^Eu, the phases of the order pa- 
rameter on superconducting islands (measured with re- 
spect to the lower superconductor) are pinned to the 
value 4>k = 0. We distinguish two different operating 
regimes of the device: one for the braiding procedure 
and one for initialization and readout. 

Flux-controlled braiding. During the braiding pro- 
cedure we set $0 = so that the charging energy of the 
large island can be completely neglected. The charging 
energies of the small islands can be considered perturba- 
tively [17], resulting in long-range Coulomb couplings. 
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between the Majorana bound states in the corresponding 
island [21]. The offset charge qk accounts for the effect 
of nearby gate electrodes. If Uk ^ Em, the low-energy 
sector of the system is described by the effective Hamil- 
tonian (see Appendix A) 

-ffbraiding = -iAijBlE ^ i/^2lElF " iAsjElC, (2) 
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where ai2 = (e/2;i)($i-h$2), "23 = (e/2;i)(<I>2-l-$3), and 
031 = —0L\2 — a23 are gauge-invariant phase differences 



between the smaller islands. The three couplings A^ are 
all tunable with exponential sensitivity via the fluxes $i, 
increasing from A^i^ (the off state) to Amax (the on 
state) when increases from to <I>inax < h/Ae. On 
the other hand, the tunnel couplings at the T-junction 
vary slowly with the fluxes, so the three overlapping Ma- 
joranas remain strongly coupled throughout the opera- 
tion. 

Out of the six Majorana operators, we define three 
fermionic creation operators c\ = -^{ja + t-Jb), = 

5(70 + «7d), and c\ = ^{'^e + ilp)- We wfll braid 
the Majoranas 73 and 7c by using 7^ and 7^ as an- 
cillas, as specified in Fig. 2. At the beginning and at 
the end, the Majoranas 75 and 7f are strongly cou- 
pled (|$2| = •i'max)- If all other couplings are off we 
are left with two degenerate states that define a topolog- 
ical qubit. In the odd-parity sector they are (p) = |10)|0) 
and (°) — |01)|0). During the exchange of Majoranas 75 
and 7c the fluxes $1, $2, ^1^3 are varied between and 
i'&inax according to the table shown in Fig. 2b. Com- 
puting the non-Abelian Berry phase for this adiabatic 
cycle as in Ref. [22] shows that braiding has the effect of 
multiplying the topological qubit state with the matrix 
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up to corrections of order Amin / Amax , with 
/A„iax < 1 because of the exponential sensi- 
tivity of these quantities on magnetic fluxes. Repeating 
the cycle n times corresponds to applying the gate U". 

Initialization and readout. Braiding needs to be 
preceded and followed by a readout of the topological 
qubit. For that purpose, before and after the braiding 
flux cycle we increase $0 from to <l>max, so that the 
spectrum of the transmon depends on the fermion parity 
V — ilAlB [21]. After having initialized the ancillas in 
the state ]0), we set $1 = $2 = 'i's = 0, to decouple 
the four Majoranas 7c, 7r>, 7b, 7f from 7^,75 and to 
minimize the effect of residual cross-capacitances [32]. 

In this configuration it is possible to execute a projec- 
tive measurement on the fermion parity V by irradiating 
the resonator with microwaves. The system composed 
by the transmon qubit and microwave resonator can be 
described by the Hamiltonian 
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Here, uq is the bare resonance frequency of the cavity, g 
is the strength of the coupling between photons and the 
transmon qubit, and MIq ~ y^SEj^Ec is the transmon 
plasma frequency, with Eq the charging energy of the 
transmon including the contributions of the small islands. 
We have defined a± = {a.^ i ia-y)/2 and 
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Figure 3. Panel a): Minimal transmon circuit for fully flux-controlled topological qubit. The nanowires are placed in a 
triangular loop formed out of three T-junctions [28]. In this geometry, all single-qubit Clifford gates can be implemented. 
Panel b): Schematic overview of a Random Access Majorana Memory consisting of eight topological qubits. Compensating 
fluxes (dotted circles) are included between the topological qubits to ensure that the gauge-invariant phase differences in the 
different topological qubits are independent of each other (see Appendix C). 



where 6ei, Seq cx exp{~ ^/SEJ^JEc) are determined by 
the energy levels £„ = £"„ — (— 1)"&„ cos(7r(7o/e) of the 
transmon [17]. We assume that the induced charge is 
fixed at go = for maximal sensitivity. 

The transmission line resonator is typically operated 
far from resonance, in the so-called dispersive regime 
[17, 19, 20], when {n + l)g^ <^ Suj'^, with n the number 
of photons in the cavity and Slu ^ D,q — ujq. The Hamil- 
tonian (5) then produces a parity-dependent resonance 
frequency (see Appendix B) 
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A flip of the topological qubit can thus be measured as 
a shift in the resonance frequency by the amount 
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The probability of observing a change in the resonance 
frequency of the cavity after n consecutive braidings, 
Pfiipin), is dictated by the Majorana statistics: pmp^n) — 
|(1|Z^"|0)|2 = |(0|Z^" |1)|2. The sequence of probabilities, 
Pflip = IjIjIjO foi' ^ = li2, 3,4, repeats itself periodi- 
cally. Therefore, the non-Abelian nature of Majoranas 
can be probed by collecting statistics for different values 
of n. 



II. RANDOM ACCESS MAJORANA MEMORY 

The TT-circuit of Fig. 2 is the minimal circuit which 
can demonstrate non-Abelian Majorana statistics, but it 
does not allow for the application of two independent 
braidings. The full computational power of Majoranas 
can be achieved by increasing the number of T-junctions. 
We adopt the triangular loop geometry introduced by 



Sau, Clarke, and Tewari [28], which is the minimal circuit 
for a fully flux-controlled topological qubit (see Fig. 3a). 
It consists of six Majorana islands placed between the 
upper and lower superconducting plates of a transmon 
qubit, referred to as bus and (phase) ground respectively, 
and a transmission line resonator for the readout. 

In this geometry the braiding and readout can be per- 
formed in a similar way as in the case of the 7r-circuit. In 
the braiding configuration, we set $o = 0- Any pair of 
the Majoranas 7a, 7b, 7c j and 713, forming the topolog- 
ical qubit, can now be braided with the help of magnetic 
fluxes {k = 1,2,..., 6), while and 7f serve as an- 
cillas. Moreover, the parity of any pair of Majoranas 
can be measured by moving them to the "measurement" 
island, the one coupled to the bus via the flux $1 in 
Fig. 3a. During the measurement $fc (fc = 1, 2, 6) 
and $0 = 3'max, so that all the small islands are coupled 
via large Josephson energy either to the bus or to the 
ground. Therefore, the measurement configuration is de- 
scribed by the readout Hamiltonian (5), where V is the 
parity of the two Majoranas in the measurement island. 

Since the typical length of a transmon is hundreds of 
microns, it is in principle possible to scale up the design 
by considering a register of several topological qubits, 
shown in Fig. 3b. The measurement configuration is 
still described by the readout Hamiltonian (5) (see Ap- 
pendix C), where the parity operator is now 
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Here 7„x and 7„y denote Majorana fcrmions on the mea- 
surement island belonging to topological qubit n. Thus, 
a readout of the resonance frequency corresponds to a 
projective measurement of this multi-qubit operator. Al- 
though the product in Eq. (8) runs over all N qubits, we 
can still choose not to measure a qubit by moving the 
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Figure 4. Quantum circuits for universal quantum computation in the ramm. In this figure, pi,p2,P3 ~ ±1 represent results of 
projective single- or multi-qubit measurements, whose outcomes, carried by classical channels (double lines), determine post- 
selected unitary operations. Panel a): CNOT gate. Here Ri = exp [ija^il — pi)], R2 = exp [ifp2P30"z], -R3 = exp [i^P2P30-x], 
R4 = exp [—ijPso-x] are all gates obtainable by braidings. Panelb): 7r/8 phase-gate T = diag (l, exp i^) , relying on distillation 
of the state |^) = (|0) -|-expij|l)) /V^. The required unitary operations are in this case R^ = exp [—i^az{l — pi)] and 
Ra = Ri- Panel c): teleportation protocol. Here R — exp [ijaz{l — P1P2)] exp [ija^il — ps)]. Apart from teleporting the 
unknown quantum state the protocol leaves the remaining two qubits in an entangled Bell state 



corresponding pair of coupled ancillas jnE , 7nF to the 
measurement island. Because these ancillas are always 
in a state |0), they do not influence the measurement 
outcome. Since the Majorana fermions can be selectively 
addressed, we call this architecture a Random Access Ma- 
jorana Memory (ramm). 

The number of qubits in a ramm register cannot be in- 
creased without limitations. Firstly, the frequency shift 
Wshift decreases with the number of topological qubits. 
The main decrease is caused by the reduction of the cou- 
pling with the number of topological qubits, which 
occurs because the Majorana fermions at the T-junctions 
are localized in three different islands (see Appendix C). 
An additional decrease is caused by the renormalization 
of the total capacitance of the transmon due to the small 
islands. Furthermore, each topological qubit introduces 
an extra pathway for quasiparticles to be exchanged be- 
tween the bus and the ground. Such quasiparticle poi- 
soning rates at thermal equilibrium are negligibly small 
and the poisoning due to non-equilibrium quasiparticles 
can, at least in principle, be controlled by creating quasi- 
particle traps. 

The limited number of qubits is not an obstacle for 
the scalability of quantum computation. Beyond this 
limit, the computation can be scaled up by using several 
transmons in a single transmission line resonator, and the 
coupling between the topological qubits in different reg- 
isters can be achieved by introducing tunable Josephson 
junctions between the transmons. Furthermore, the com- 
putation can be parallelized, because transmons can be 
coupled to several different transmission line resonators 
[33-35]. 



III. MULTI-QUBIT MEASUREMENTS AS A 
SOURCE OF COMPUTATIONAL POWER 

Multi-qubit measurements in the RAMM offer two sig- 
nificant benefits. Firstly, these measurements can be ap- 



plied without any locality constraint, so that the quan- 
tum fan-out [34] , the number of other qubits with which a 
given qubit can interact, can become large for the ramm 
architecture. Secondly, the overhead in the computa- 
tional resources can be reduced because the products of 
Pauli matrices involving several topological qubits can be 
measured directly. We demonstrate these advantages in 
the realization of a universal set of gates, fast creation of 
maximally entangled states, and implementation of error 
correction schemes. 

Quantum gates. All single-qubit Clifford gates, the 
CNOT gate, and the tt/S phase gate required for univer- 
sal quantum computation [36], can be realized in the 
RAMM with errors that are exponentially small in macro- 
scopic control parameters (see Appendices B and D). 
Single-qubit Clifford gates can be realized with braid- 
ing operations only, and the quantum circuits for the 
two remaining gates are summarized in Fig. 4. The 
CNOT gate, shown in Fig. 4a, is a modified version of the 
Bravyi-Kitaev algorithm [37, 38] involving three topologi- 
cal qubits (target, control, and one ancilla). Efficient tt/8 
phase gate implementations are based on distillation pro- 
tocols [39], requiring several noisy qubits to prepare one 
qubit in a particular state \A) = (|0) +e*''/'*|l)) /V2. 
This state can then be used to perform the tt/8 gate using 
the circuit shown in Fig. 4 b. Distillation may take place 
in dedicated ramm registers (see Appendix D) in par- 
allel with other computation processes, and the distilled 
state can be teleported to the computational register (see 
Fig. 4c). 

Preparation of 2D cluster states. The ramm can 
be used to efficiently create maximally entangled multi- 
qubit states, such as 2D cluster states [40-42], which 
make it possible to realize any quantum circuit by means 
of single-qubit operations and measurements [43]. 

To generate a 2D cluster state in the RAMM architec- 
ture one has first to assign a label to each topological 
qubit in order to establish its position and neighbours 
on a logical lattice (see Fig. 5a). Due to the non- locality 
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Figure 5. Preparation of a 9-qubit 2D cluster state with a 
RAMM. The nine qubits (represented by circles) are arranged 
in a 3 X 3 square logical lattice, and numbered from left to 
right and top to bottom. Panel a): The nine stabilizer op- 
erators Ki, . . . Kg necessary to prepare the 2D cluster state. 
They are products of Pauli matrices, involving all qubits con- 
nected by lines, with black and grey dots representing ax and 
(Tz operators, respectively. Panel b): The quantum circuit 
creating the 2D cluster state in a 9-qubit ramm register, con- 
sisting in a sequence of projective multi-qubit measurement 
of the 9 stabilizers. 



of measurements in the ramm, the logical lattice does 
not need to be related to the physical system. The clus- 
ter state may be prepared in several ways [40, 42]. An 
efficient procedure requires measuring the stabilizers 



(9) 



where a goes through all sites of the logical lattice and (i 
labels the nearest neighbours of a. The total number of 
measurements required is equal to the number of qubits 
in the cluster state. In Fig. 5b we draw a circuit to create 
the 9-qubit 2D cluster state in a RAMM register. To verify 
their entanglement properties, one possibility is provided 
by the teleportation protocol of Ref. 43. 

Efficient quantum error correction. Although 



topological qubits have intrinsically low error rates, 
grouping them into a ramm register allows to addition- 
ally implement efficient error correction. Error correc- 
tion schemes [36, 44] are based on measurements of stabi- 
lizer generators, which are products of Pauli matrices be- 
longing to different qubits. The measurement outcomes 
give error syndromes, which uniquely characterize the er- 
rors and the qubits where they occurred. The ramm 
allows for efficient error correction schemes, due to the 
possibility of measuring stabilizers of different length, as 
well as correcting errors using single-qubit Clifford gates. 
There are two advantages in comparison with architec- 
tures where only single- and two-qubit operations are 
available: higher error thresholds and reduced overhead 
in computational resources. 

In order to quantitatively compare these advantages, 
we consider the 7-qubit Steane code [45] as a concrete 
example of quantum codes, and assume a realistic error 
model. We ffiid that the error threshold of the ramm can 
be an order of magnitude larger than the error threshold 
of a reference architecture that can only perform single- 
and two-qubit operations (see Appendix E). Additionally, 
the ramm implementation of the Steane code is much 
more compact. Already in the first level of concatenation, 
the fault-tolerant implementation of syndrome measure- 
ments in the reference architecture requires 24 ancillas for 
each logical qubit, while none are needed in the ramm. 

Although we have calculated the improvements only 
for the 7-qubit Steane code, the advantages are charac- 
teristic for all error correction schemes, including surface 
codes [46, 47]. 



IV. DISCUSSION 

To control and manipulate quantum information con- 
tained in the Majorana zero-modes of superconducting 
nanowires it is necessary to braid them and measure 
their parity. We have designed a transmon circuit where 
both operations can be performed by controlling the 
magnetic fluxes through split Josephson junctions, with- 
out local adjustment of microscopic parameters of the 
nanowires. The minimal circuit for the demonstration 
of non-Abelian Majorana statistics is a 7r-shaped circuit 
involving four independent flux variables. An extended 
circuit consisting of many topological qubits in parallel 
allows for non-local multi-qubit measurements in a Ran- 
dom Access Majorana Memory, providing the possibil- 
ities of efficient creation of highly entangled states and 
simplified (ancilla-free) quantum error correction. 

Since all the requirements for the realization of the tt- 
circuit and ramm are satisfied with the typical energy 
scales of existing transmon circuits and transmission line 
resonators (see Appendix F), flux-controlled circuits are 
a favorable architecture for the demonstration of non- 
Abelian Majorana statistics and the realization of fault- 
tolerant quantum computation. 
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Figure 6. The 7r-shaped transmon circuit discussed in the 
main text, reproduced here with labels of the ten Majorana 
bound states. 
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Appendix A: Theoretical description of the vr-shaped 
circuit 

The TT-shaped circuit discussed in the main text is re- 
produced here in Fig. 6. We label the two superconduct- 
ing plates forming the transmon "bus" and "ground" , 
both hosting two Majorana bound states, labeled 7bi, 7f,2 
and 7gi,7g2 respectively. The smaller superconducting 
islands are labeled with an integer k = 1,2,3. Each of 
them supports two Majorana bound states 7fci,7fc2- We 
will work in a gauge where all phases are measured with 
respect to the phase of the ground island. We denote 
with (p the phase of the bus and with (/j^ that of the k-th 
island. 

We start from the Lagrangian of the system, 



£ = T-Vj~Vm. 
The first term is the charging energy 
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k = l 



(A2) 



Here Cq is the capacitance between bus and ground, while 
Cg,/c {Cs.k) is the capacitance between the fc-th Ma- 
jorana island and the ground (the bus). The last two 
terms include the induced charge go on the bus and qk 
on Majorana islands. The effect of cross-capacitances be- 
tween Majorana islands is negligible assuming that they 



are small in comparison with the capacitances to the bus 
and the ground. 

The second term is the Josephson potential 



Vj^Ej^Q{^){l-cosi 



3 

fc=i 



j,fc($fe)(l - cos (t)k). 



The Josephson energies £'j,o(*i'o) = 2£^j_o(0) cos(e<i>o//i) 
and Ej^ki^k) — 2i?xfc(0) cos(e<I>fe/fi) can be varied 
in magnitude by changing the fluxes between and 
I'i'maxI ^ h/Ae. We are assuming for simpHcity that the 
split junctions are symmetrical, but this requirement can 
be removed without affecting our results. 

The third term is the Majorana- Josephson potential 



Vm = Em [ilb2lgi cos (i0 -I- a^^g) 
+ ilgilii cos (Ofgi - 
+ j7ii7b2 cos + aib)] 

+Em [1712721 cos (5 01 - + ^12) 
+ 1721731 COS [\4>2 - \cj)3, + 0123) 

+ 1731712 COS (i^a - \cj)i + azi)] ■ 



(A3) 



The two square brackets in this expression group the 
terms corresponding to the two T-junctions. All tun- 
nel couplings are for simplicity assumed to be of equal 
strength Em- The arguments of the cosines include 
single-electron Aharonov-Bohm phase shifts between dif- 
ferent islands, 

abg^e<^Q/2h (A4a) 

agi = e<^i/2h (A4b) 

aib = - (e$o + e$i)/2;i (A4c) 

"12 = (e*i + e<i>2)/2h (A4d) 

"23 = (e*2 + e$3)/2/i (A4e) 

aai = - (e$i -I- 2e<I>2 + e$3) 72/1 (A4f) 

There is a constraint between the charge contained in 
each superconducting island and the parity of the Ma- 
jorana fermions belonging to that island [48]. The con- 
straint can be eliminated via a gauge transformation [49] 



k=l 



^ink<t>k/2 



J«7bl7ft2 , nk 



j«7fci7fc2 • 



(A5) 
(A6) 



where the products extends over all Majorana junctions. 
The transformation has two effects on the Lagrangian: 

• it changes the induced charges appearing in Eq. 
(A2), 



90 



qo + en , qk -> qk + euk 



(A7) 



so that the Majorana operators enter explicitly in 
the charging energy, and 
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• it modifies the Majorana-Josephson potential 
fl^VM^^ so that it becomes 27r-periodic in all its 
arguments 4>k ■ 

In the following, we will work in this new gauge where Eq. 
(A7) holds. The explicit form of VI^Vm^ is not necessary 
here, as we will only need the equality 



Vm\ 



<Pk=4>=Q 



(A8) 



which is trivial since £71 



= 1. 



10, =0=0 - - Starting from 
the Lagrangian (Al), we will now derive the low-energy 
Hamiltonians used in the main text for the braiding and 
the readout. 



1. Braiding 

When we want to braid or move the Majoranas, we 
maximize the energy i?j.o('I'o) by setting $o = and we 
require the condition 



Ejfi{Q),Ej^k{'^k) » Em,Ec,Ecm 



(A9) 



where Ec,o = IIC^ and Ec,k = /'2{CB,k + Cgm)- 
Since the Josephson term Vj dominates over the kinetic 
and Majorana terms T and Vm, the action S = J f-dt 
is then minimized for (f) ~ (j)^. — and (f> — 4>k ^ 0- 
All the superconducting islands are in phase. Under the 
additional condition 



Ej,o{0) ^ Ej^ki'i'k) 



Ec.o 



Ec,i 



(AlO) 



we can neglect quantum phase slips around the minimum 
= 0, but not around the other minima (pk = 0. The low- 
energy Hamiltonian Hyi then contains only the Majorana 
operators: 



3 

E 

fc=i 



(All) 



In this passage we have introduced the six Majorana op- 
erators "fAnB,lCilD,lE,lF, given by 

7A=7fei, (A14a) 
_ cos agi7fc2 + cos aiblgi + cos abglii (^p^i^^,-) 
cos^ Ugi + cos^ aib + cos^ abg 

(Al4c) 
(Al4d) 

(A14e) 
(A14f) 



7c = 732, 
ID = 7s2, 



IE 



cosa237i2 + cosa3i72i -I- cosai273i 



VcOS^ q;23 -t- COS^ 031 -|- COS^ ai2 
7-F = 722- 

The coupling strengths are 

cos ajjg 



Al = C/i 



yjcos^ agi + cos^ aib + cos^ abg 
cos a23 

-y/cos^ ai2 + COS^ a23 + cos^ a3i 
COS 0131 



' Vcos^ ai2 + cos^ a23 -I- cos^ aai 
A3 = 

VCOS^ Q!l2 + COS^ a23 + cos^ a3i 



(Al5a) 
(Al5b) 
(Al5c) 



2. Readout 

During the readout of the transmon qubit, we set $o = 
$max, so that the Josephson energy Ej^ is minimized, 
and all = 0. We require then that 



Ej,k{0) £^J,0($max) 

Enk 



Ec.o 



(A16) 



In physical terms, all Majorana islands are now in phase 
with the ground: (j)k = (pk = 0. Neglecting quantum 
fluctuations and phase slips around these minima, we 
may re- write the Lagrangian in a form that depends only 
on (h 



where 



Uk = 16 



Ec,kEl^^ 



27r2 



^-yJ&E.^JEck cos(gfc7r/e). 



(A12) 

is the tunneling amplitude of a phase slip process from 
(/>A: = to = ±27r [17], also reported in Eq. (1) of the 
main text. 

There are still ten Majorana operators in the Hamil- 
tonian (All), but we can eliminate four of them by as- 
suming that the tunnel couplings are stronger than the 
Coulomb couplings: Eyi 3> Uk- To first order in pertur- 
bation theory in the ratio Uk/Eyi, we then obtain the 
Hamiltonian used in the main text 



Ejfl (1 - cose/)) - VI^VkM 



t=0 



(A17) 



Apart from the contribution of the term Vm, the whole 
system can be treated as a single hybrid top-transmon 
[21], with Josephson energy Ej^ and capacitance 



C — Co + ^ CB,k- 



(A18) 



fe=i 



In the regime Ej^q ^ Ec — 6^/2(7, the energy levels of 
the transmon are given by [17] 



H = ~iAi7B7_E - i/^2lElF - i^slElc 



(A13) 



-l)"(Se„i76i7b2COs(7rg/e) , 



(A19) 
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where 



are, respectively: 



(A20) 



21 I 3 



^/SEj,„/Ec ^ (^21) 



Taking into account the two lowest levels of the trans- 
mon (n = 0, 1), we arrive at a low-energy Haniiltonian 

^^top-transmon = Cr^ [11^0 + i-fbllb2 5+ COs(7rgo/e)] 

+ilbilb2 5- cos(7rgo/e) + ^^Vm^^^^^^^^ (A22) 

with definitions MIq — £i ~ ep: 5± — (fci ± fco)/2. The 
Pauli matrix acts on the qubit degree of freedom of 
the transmon. For 5± <C Em, the low energy sector of 
this Hamiltonian can be written in terms of 7a, ■ ■ ■ ,7_f 
as 

-fftop-transmon = <J z [j^^^O + ^lAlB A+ C0s(7rqo/e)] 

+ ilAlB A- cos(7rgo/e) (A23) 



where 



Ah 



5± cos Ugl 



cos^ abg + cos^ agi + cos^ 



(A24) 



When combined with the Jaynes-Cummings Hamiltonian 
describing the coupling with the resonator, this Hamil- 
tonian reproduces Eq. (5) of the main text. The inter- 
action with the microwaves will be described in detail in 
the next Appendix B. 



Appendix B: Measurement through photon 
transmission 

The Hamiltonian i?rcadout of the main text describes 
the coupling between the top-transmon and the cavity 
modes in the system through a Jaynes-Cummings inter- 
action of strenght g. In particular the fermionic parity 
of the transmon is a conserved quantity in the Hamil- 
tonian whose energy levels will directly depend on the 
value of V. 

We assume that the induced charge is fixed at go = 
to maximize the sensitivity of the read-out. The 
Jaynes-Cummings interaction couples the pairs of states 
(|n,t, "P), \n + l,i,V)) where n and n + 1 label the num- 
ber of photons in the cavity and | t) , | i) denote the two 
lowest energy eigenstates of the transmon. Therefore, 
the eigenstates of ^^roadout are in general superpositions 
of the kind ajn, f, "P) +/3|?^+l, i, "P) with the exception of 
the uncoupled vacuum states |0,4,, P). Their eigenvalues 



£n,±,v = + ^ j ^'^o + PA_ + 

± ]^^/{h5oj + 2PA+)2 + Ah^g"^ (n + 1) (Bl) 



eo,p -P(A_ - A+) - -nrtn. 



(B2) 



In the dispersive regime, Slo'^ ^ g'^{n + 1), the energies 
£n.±.v can be approximated at the first order in g'^/Suj'^ 
as: 

. ^ h^g^(n + l) 

= nnu:o + P (A_ + A+) + -H^o + jr^f;^^ 

(B3) 

en+i,i,v = (n + l)hojo + P (A_ - A+) + 

1 h^g^ (n + 1) 



2 h6io + 2PA I 



(B4) 



The respective eigenstates are approximately 171,1, P) 
and \n + 1,1, P) up to corrections of the order g'^/Suj'^. 
From the previous equations it is easy to obtain the ef- 
fective resonance frequency lUcS (P) and its shift Wshift 
corresponding to the different states of the topological 
qubit. Since we are considering the dispersive regime 
with a positive detuning, ilg > ujq, we assume in the 
following that the state of the transmon remains in the 
ground state | |). 

We also point out that in the Hamiltonian -ffrcadout we 
are neglecting the excited states of the transmon, which 
result in a renormalization of the parameters, including 
<^shift: through virtual transitions. The precise expres- 
sions for the renormalized parameters are known [17], 
but are not needed here. 

To perform the measurement of the topological qubit 
we introduce in the cavity photons with a frequency 
which is approximately Wcff(P = +1)- The photon trans- 
mission probability T_|_ for the state |P — 1) is then larger 
than the probability T- corresponding to |P = —1). We 
count the number of photons riph that passes through the 
cavity during a measurement time t^- The probability 
distributions for nph in each state are Poissonian, and for 
sufficiently long measurement time can be approximated 
with normal distributions 



P(nph, |P = ±1)) = Pois(nph, A±) « N(nph, A±, y/X±) 

(B5) 

where X± oc T±tMH and k ~ 1 — 10 MHz is the cavity 
decay rate. Since T+ > T_, also A+ > A_. 

We decide that the measurement outcome is P = +1 if 
riph > a; = ^/\+\Z and the outcome is P = — 1 if riph < 
X. Therefore the error of the measurement outcome is 
given by the following; 



1 



dn ( —in — A I )2 

exp ' 



dn 



2 J, ^2^ 



: exp 



2A+ 

2A_ 



(B6) 
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Since A + , A_ 1 



fermions 



where 



2a; -^/tt' 



/A7 - ,A 



(B7) 



V2 



(B8) 



We notice that the probability of a measurement error 
decreases exponentially with utM- On the other hand, 
the probability of storage error, namely the chance that 
the topological qubit will decay during a time interval ^a/, 
increases as A^intAi /fi. Because Amin/K can be made 
exponentially small in macroscopic control parameters, 
exponentially small measurement errors can be achieved. 



Appendix C: Low energy Hamiltonian for a Random 
Access Majorana Memory architecture 

We will now describe an effective Hamiltonian for 
RAMM architecture hosting N topological qubits, such as 
the one shown in Fig. 3 of the main text. Fig. 7a shows an 
equivalent setup, including only two topological qubits. 
By including compensating fluxes 



comp,n 



(Cl) 



fc=i 



after each topological qubit, the gauge invariant phases 
in each topological qubit are independent of each other. 
The single-electron Aharonov-Bohm phase-shifts an,kk' 
at the tunnel junction between islands k and k' of the 
n-th qubit are then given by 



a„,i2 = e($o + + *n,2)/2/i 

a„,26 = e($„,2 + 2$„,3 + 2$„,4 + 2$„,5 + ^n,6)/2h 
anfii = -e($o + + 2$„,2 + 2$„^3 

+ 2$„,4 + 2$„,5 + $„,6)/2S 

a«,23 = e($„,2 + $n,3)/2fi 
a„,34 = e($„,3 + $„,4)/2;i 

a„,42 = -e{<^n.2 + 2$«,3 + *«,4)/2/l 

a«,45 = e(<I>„,4 + $„,5)/2;i 

0^71,56 



= e($„,5 + $„,6)/2/i 

= -e($„,4 + 2$„^5 + *«,6)/2/i. 



(C2) 



By starting from a Lagrangian and following a simi- 
lar approach to that of Appendix A, we find that the 
low-energy Hamiltonian is described by six Majorana 



ln,A — 7n,32, 

COSa„,347„,22 + cos Q;„,427n,31 + COS an,237n,41 

In.B = 



ln.,C 



^/cos2 a„,23 + cos2 a„_34 -I- cos2 a„_42 

COS an^567n,42 + COS Q;n,647n,51 + COS an,457n,62 
^COS^ an,A5 + COS^ Q!„,56 -|- COS^ Q!„,64 



7n,-D — 7n.52, 

COSa„,267n,12 + COS a„.6l7Ti,21 + COSa„,i27n,61 

In.E = 



^ cos2 a„,i2 -t- cos2 a„,26 + cos2 a„,61 

7n,F = 7n,ll (C3) 

that form the triangular loop network of Fig. 7b. 



1. Low-energy Hamiltonian in braiding 
configuration 

In the braiding configuration (f>o — Oi and the low- 
energy Hamiltonian is, for each qubit rt. 



(n) 
qubit 



il\n,\lFlE - iAn,2lElB - iAn.slBiA 

(C4) 



where long-range Coulomb couplings A„jt are 
/V _jr cosa„,26 

A,i^2 = Un,2 



An,3 
A„^4 



^/cos^ a„42 + cos^ an,26 + cos2 a„ 61 

COS a„,34 

Y^COS^ a„,23 + C0S2 a„,34 -I- C0S2 a„,42 

COS a„,6i 

^COS^ a„^i2 + COS^ an,26 + C0S2 anfil 
COS a„,42 



^/C0S2 a„ 23 + C0S2 a„ 34 -I- COS^ a„,42 

COSQn^23 

^/cOs'^ a„,23 + COS^ a„,34 -I- C0S2 a„,42 
COS an, 56 



in, 6 



Un,5 
Un.6 



y'cOS^ anA5 + COS^ a„,56 -t- COS^ a„,64 
COS Q!„_46 



^/cos2 a„,45 -I- cos2 a„,56 -I- cos2 a„,64 

COS an,12 

Y^cos^ a„42 + cos^ a„,26 + cos^ a„,6i 

cos 0;n,45 

^C0S2 a„_45 -f C0S2 a„_56 -I- COS^ Q!„_64 



(C5) 



For computational purposes, one should be careful that 
the An,k do not change signs during the variation of 
the magnetic fluxes that takes place during a compu- 
tational process. This may happen if some of the an,kk' 
in Eq. (C2) cross the value 7r/2. However, during any 
computation, maximally two of the fluxes are simultane- 
ously turned on. Therefore, it is always possible to adapt 
the signs of the magnetic fluxes in such a way that the 
fluxes can be tuned in a range |$n,A:| = [0, $max], where 
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b) 



7FO 
Ai- 

A2/ 



IB i 

A3 



A4 



.A« 



»7c 

I As 



Figure 7. Panel (a): Part of the ramm circuit showing two fully-controllable topological qubits. Compensating fluxes are 
included between the topological qubits in order that the gauge-invariant phase differences in the different topological qubits 
are independent of each other. Panel (b): opological qubit formed by the six Majorana fermions. The six couplings Ai, . . . , As, 
see Eq. (C5), can all be individually controlled by the fluxes $1, . . . , $6. The parity of the two Majoranas coupled by Ai can 
be measured, as explained in Appendix C 2. 



'I'max < h/Ae. We also notice that the signs of the cou- 
phngs Aj in Eq. (C4) depend on the signs of the micro- 
scopic tunnel couplings Em- These signs will determine 
the chirality of the braiding of the Majorana fermions in 
each T-junction. 



2. Low-energy Hamiltonian in the readout 
configuration 

During the readout, we set $0 = ^'max and all other 
fluxes ^n,k = 0. Following the same reasoning of Ap- 
pendix A 2, we set (j)n,i = 4> and 4'n.,k=^i = for each topo- 
logical qubit. The Lagrangian for the RAMM becomes 



2e 



N 



0„,fc=O 



(C6) 



where vj^^'' describes the Majorana- Josephson potential 
for the three T-junctions in each topological qubit n, 

6 



= n e'^'" 



■i7n,fcl7Ti,fc2)0fe/4 



(C7) 



fc=i 



N 6 



N 



C^Co + J2J2 ^B,k + J2 (C8) 



n=l k=2 



and 



N 



9tot = 90 + ^ qn,l- 



(C9) 



The low-energy Hamiltonian of the system can now be 
derived analogously as in Appendix A 2. By using the 



equality 



N 



COS |^7rgtot/e + tt ^ (5-5 «7n,ii7n,i2) J 

N 

= W «7n,ii7n,i2 cos(7rgtot/e) (CIO) 

n=l 

we find 

^RAMM = cTz [\ti^o + A+ cos(7rgtot/e)] 

+ 7^ A_ cos(7rgtot/e) (CU) 

where V is now the joint parity operator of the Majorana 
fermions at the measurement islands 



N 



"P = iln.Fln,E- 



(C12) 



n=l 



The couplings A± decrease exponentially with the 
number of topological qubits involved in a single RAMM 
register 



A. 



N 

n 



cos a„ 96 



1 



/cos^ a„,i2 + cos^ a„^26 + cos^ a„^6i 

(C13) 

In the design of a ramm register, shown in Fig. 3b in 
the main text, the frequency shift Wshift is decreased by 
all topological qubits, including the ones which are not 
involved in a given multiqubit measurement. This limi- 
tation of RAMM can be relaxed in a more optimal design, 
where additional tunable Josephson junctions are intro- 
duced from the measurement island to the ground. In 
this case only the topological qubits, which are involved 
in the given measurement, contribute to the decrease of 
frequency shift. The expense one needs to pay for in- 
troducing new Josephson junctions is that the gauge in- 
variant fluxes have more complicated magnetic flux de- 
pendence and several Josephson couplings need be si- 
multaneously controlled when the Coulomb couplings are 
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turned on. We point out that although we have exphcitly 
considered the control of the Coulomb couplings with the 
help of magnetic fluxes, at least some of the macroscopic 
control parameters Ej^k/Ec,k of the superconducting is- 
lands can alternatively be controlled with gates. 



Appendix D: Universal gates for quantum 
computation 

The RAMM setup allows us to perform universal quan- 
tum computation in a fault-tolerant way. To show this, 
it is necessary to implement a universal basis of quan- 
tum gates using only braiding operators and multi-qubit 
measurements as building blocks, thus ensuring the pos- 
sibility of obtaining arbitrary multi-qubit gates with er- 
rors that are exponentially small in the macroscopically 
tunable parameters. One possible set of gates allowing 
for universal quantum computation are the single-qubit 
Clifford gates, the CNOT gate and the 7r/8 phase gate. In 
the following we explain how to realize these gates in a 
RAMM architecture. 



1. Notation 

Each topological qubit n has four computational Majo- 
ranas 7n,A) In.S) 7n,Cj 7ri,D Eind two ancillary Majoranas 
ln,E,ln,F, which are needed to move or braid the com- 
putational ones. The Pauli matrices for each qubit can 
be chosen as 

<yn,z = iln.Aln^B (Dla) 
(yn,x = iln,Bln,C (Dlb) 
0-„,y = iln.Aln^C- (I^lc) 

2. Single-qubit operations 

Projective measurements on the Pauli basis and a set 
of Clifford gates can be obtained by manipulating the 
positions of the four computational Majorana fermions 
in the triangular loop geometry. The positions of the 
computational Majoranas 7n,A, 7n,B, 7n,Ci 7n,_D can be 
changed using the ancillary Majorana jn,E,"fn.F, which 
remain strongly coupled throughout the process. The 
corresponding qubit transformation can be derived either 
by a direct computation of the non-Abelian Berry phase 
acquired by the ground state wave function of the Hamil- 
tonian (C4), or by following the evolution of the Majo- 
rana operators in the Heisenberg picture, as explained in 
detail in Ref. [31, 50]. 

Exchanging the positions of 7„,yi,7n,B (as represented 
in Fig. 8a) or jn,B and jri,c (Fig. 8b) respectively yields 



the braiding gates 

= e-^f , (D2) 
= e^^?'^- . (D3) 

The chirality of the braiding operations (i.e., the sign 
of the exponent in lAzMx) is determined by the signs 
of the couplings of the qubit Hamiltonian, Eq. (C4). 
Physically, the sign depends on the induced charges on 
the Majorana islands, the values of the fluxes and the 
signs of the microscopic tunnel couplings ±i5jvf at the 
T-junctions. Here, we have made a specific choice of 
chiralities. Another possibility of chiralities would not 
be harmful as long as they remain constant during the 
computation processes. 

A combination of these two operations yields the quan- 
tum gate corresponding to the braiding of and 7c 

Uy ^ Ul ^ e'il'^y . (D4) 

When combined with the 7r/8 phase gate described in Ap- 
pendix D 4, these quantum gates are sufficient to realize 
any single-qubit rotation. 

To realize projective measurements on an,z (or <Jn.x)^ 
we first need to bring the two Majorana fermions 
ln,A,ln,B (or 7n,B,7n,c) ou the island connected to the 
bus, the one occupied by ^n.E,ln,F in Fig- 7a. Then we 
measure the fermion parity operator (C12), where now 
the two Majoranas ^n.E,ln,F are replaced by the compu- 
tational ones. For instance, in the case of a measurement 
of cr„^2, we would measure the operator 

V = i^n,Aln,BW_ilk.Elk,F = Crn,z, (D5) 

since the parity of the ancillary Majorana of each topo- 
logical qubit is preserved, Vk.EF — ilk,Elk,F = +1- In 
the end, we bring the two computational Majoranas back 
to their original place. The whole operation, represented 
in Fig. 8c and Fig. 8d for an,z and an,x respectively, cor- 
responds to the application of the projectors 

n.,„(p) = i (l+pa„,.), (D6a) 
nx,«(p) - i (1+P<T„,.) (D6b) 

to the wave function of the N topological qubits. Here, 
p = ±1 is the outcome of the measurement. Finally, a 
projective measurement on a^^y is obtained as 

nj,,„(p) = i (1 + pCTn,y) = ul Tl,,,{p)U^ . (D7) 

Multi-qubit measurements on the Pauli basis are a 
straightforward extension of these projective measure- 
ments where Majorana modes on different topological 
qubits are moved according to Fig. 8 to achieve the re- 
quired basis. 
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c) Projective measurement of o-^ 




Figure 8. Flux-controlled sequences of operations that realize single-qubit Clifford gates and projective measurement on the 
Pauli basis. 



3. CNOT gate 

Bravyi and Kitaev have demonstrated how to reahze 
the CNOT gate with an algorithm that is based on the 
following expansion [37, 38]: 



exp 



(*^7o7i7273) IV') = 26'" exp (^(1 - PiP2)7o7i) exp (^(1-^1^2)7273) exp (^-^^27275 

x^(l +P2i7274)^(l -Pi7o7i7374)IV'), (D8) 



where is an unimportant overall phase, 7^ (i = 0, 5) 
are Majorana operators and pi — ±1 are measure- 
ments outcomes. The Majoranas 74 and 75 are used 
as ancillas and the wave function is initialized in state 



(74 + «75)|V') = 0- Importantly, the Bravyi-Kitaev CNOT 
algorithm is based only on measurements and braidings 
of Majorana fermions. However, as one can see from 
Eq. (D8), its implementation requires a pair of ancil- 
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lary Majoranas that must be coupled to two computa- 
tional Majoranas in the target qubit, but must initially 
be completely independent on them. Due to the parity 
constraint in each topological qubit, this is impossible in 
the RAMM setup unless we extend the qubit layout shown 
in Fig. 3a in the main text. Rather than modifying the 
RAMM setup to account for these new ancillas, we propose 
an alternative version of the CNOT gate, which involves 



three topological qubits. This alternative version of the 
CNOT gate can be implemented with the quantum circuit 
shown in Fig. 4a in the main text. 

In this circuit the role of the first measurement, with 
result pi, and of the gate Ri is to initialize the third 
ancillary qubit in the state \0)a- After that, a CNOT gate 
with qi as a control and q2 as a target gate is obtained 
as: 



_ e^fP2P3-i,. e*fP2P3-T.,. g-iiP3.x„,. (1 ^ p^^^^^) (J ^ P2fTi,,a2,,a,.,) 92, 0) = e'' \q,,q, ffi ^2, 0). (D9) 
In terms of Majorana operators, this way of representing the CNOT relies on the following equality 

exp Q7iA7iB72B72C73yi73B) |V')i2|0)a = 2e'*exp (-^(1 + P2P3)7ia7ib) exp (-"^(1 + P2P3h2Bl2cj 

X exp (%373S73c) x(l + »P373A73C);^(1 - «P27lA7lB72B72C73S73c)|'0)l2|O)a, (DIO) 



which can be considered an extension of Kitaev and 
Bravyi result. In this case the applied projections are 
all on products of parity operators from different qubits, 
which can be reduced to the form (C12) as explained 
above (see Fig. 8); all the other operators are braiding 
operators within single topological qubits. 

4. 7r/8 Phase Gate 

To complete the set of universal single-qubit gates we 
must implement the 7r/8 phase gate 

^=Ge?f)' 

with an accuracy comparable to the other gates. 

For this purpose the best techniques are based on dis- 
tillation protocols [39] . The basic idea of the distillation 
procedure is the use of several noisy qubits to prepare one 
qubit in a particular state, \A) = (|0) + e'""/^ |1)) /%/2. A 
single ancilla qubit prepared in the state \A) is enough to 
implement the tt/S gate using the circuit shown in Fig. 4b 
in the main text. 

The distillation protocol of Ref. [39] for the state \A) 
assumes that it is possible to prepare several noisy copies 
of \A) with an average initial error < 0.14. In the 
RAMM setup this can be achieved by coupling the Majo- 
rana fermions to break the ground state degeneracy [21]. 
A single distillation step is performed starting from 15 
noisy qubits. Neglecting the errors in all the Clifford 
gates and measurements of the distillation process, the 
error of the final state after one iteration is approximately 

edist « 35e3 (D12) 

in the limit of small e^. 



Since 14 stabilizer multi-qubit measurements and 15 
CNOT gates are involved in the distillation-decoding pro- 
cedure, the error in the tt/8 gate is approximately an 
order of magnitude larger than the errors occurring in 
braiding or in a single multi-qubit measurement. More- 
over, assuming an achievable initial error = 0.01 [21] 
only a single distillation step involving 15 noisy ancillas is 
needed to achieve a final error of the same order of mea- 
surement and gate errors, estimated as Amin/Amax ^ 
10~^. If the initial errors are larger or the gate errors 
are smaller, more distillation steps and a larger number 
of ancillas are preferable. Given the amount of qubits 
required, it is realistic to imagine that the distillation 
procedure will take place in one (or several) dedicated 
RAMM registers, so that it can happen in parallel with 
all other computation processes. In this way, whenever a 
7r/8 phase gate is needed in the computation, it will only 
be necessary to teleport the distilled state |^) from the 
distillation register to the computational one. 

We also note that alternatively to the n/8 gate, the 
universality can also be obtained with the help of 7r/12 
gate. This gate can be distilled with fewer noisy copies 
of the relevant state and a single distillation step also re- 
quires less multi-qubit measurements [39] . Moreover, the 
distillation can be improved by exploiting more efficient 
error correction codes: for example in Ref. 51 a different 
procedure is proposed that enables to obtain two distilled 
states \ A) out of 10 noisy ancillas, providing a better scal- 
ing and threshold for the initial errors. Finally we must 
mention that the distillation techniques in Ref. 39 require 
not only multi-qubit measurements and braiding gates, 
but also a non-unitary dephasing process. However, it 
was shown in Ref. 52 that the dephasing process is not 
necessarily needed for the convergence of the noisy states 
to a high-fidelity final state. 
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Figure 9. Measurement of the six generators of the Steane 
code. This circuit can be realized directly in a ramm archi- 
tecture. 



Appendix E: Computation of the error thresholds 

Multi-qubit measurements give significant advantages 
in quantum error correction, as compared to the usual 
schemes where only single- and two-qubit operations are 
available. The advantages obtained are twofold. Firstly, 
multi-qubit measurements allow to significantly increase 
error thresholds. Secondly, the overhead in computa- 
tional resources can be substantially decreased. 

Quantum error correction schemes are generally based 
on measurements of multi-qubit operators, usually re- 
ferred to as stabilizer generators gi [36]. Their outcomes 
give error syndromes, Pi, which uniquely characterize the 
errors and the qubits where they have occurred. Depend- 
ing on the error correction scheme, a different number of 



errors can be corrected. 

For simplicity, we consider the Steane 7-qubit quantum 
code [45], which encodes a logical qubit into seven phys- 
ical qubits and can recover an arbitrary error occurring 
in any of the physical qubits. Its stabilizer generators are 
gi — XiX^XqXt, g2 = X2X4XQXJ, 33 — X3X4X5XQ, 
94 = Z^Zj^Zj^Z'j, 35 = Z2Z3Z5Z7, and g^ = Z^Z-iZ^^Z^. 
An error detected on the i-th qubit can be corrected by 
implementing a Xi, Zi or XiZj gate, depending on the 
type of the error. 

In order to quantitatively compare the advantages ob- 
tained with the help of multi-qubit measurements to con- 
ventional schemes, we calculate the error threshold for a 
quantum memory. The error correction circuit consists 
of periodic syndrome measurements and recoveries, in- 
terrupted by a time-interval of N time steps. Time steps 
are defined so that a single gate (or measurement) can be 
performed within one time step. Our error model con- 
sists of storage errors, gate errors, data errors during the 
measurement and errors in the measurement outcomes. 
The corresponding error probabilities are est, Cg, Cdm, and 
Eom? respectively. All the errors are considered indepen- 
dent. In order to obtain the error threshold, we need 
to calculate the probability of failure happening during 
a single period of the error correction circuit, assuming 
that no failure has happened before that point. To keep 
the calculation tractable, we assume that two errors in 
different qubits always result in failure (independently on 
the type of errors), and that this happens also when one 
of the errors occurs during the syndrome-recovery part 
of the circuit and the other error has happened earlier 
in the circuit. Moreover, we assume that the errors oc- 
curring during the syndrome-recovery part of the circuit 
never get corrected by the same syndrome-recovery part 
of the circuit. This way we obtain that the probability 
of failure during a single period of the circuit is: 



P(failurc, iV) «Pom(2) 
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Here Pom(™) is the probability of having m errors in 



the measurement outcomes. 



is the probability of ob- 



taining single error in qubit i during syndrome measure- 
ment and recovery, and P^.n = A^Cst is the probability of 
obtaining single error in qubit i during the N time steps 
between the successive error detections and recoveries. 

To estimate the error threshold we minimize the prob- 
ability of failure per time step 



Vi = min{P(failure, N) /{N 



No)}, 



(E2) 



where A'o is the number of time steps required to perform 
the syndrome measurements and the recovery. The quan- 
tum error correction threshold is obtained by demanding 
that pf = est ■ Because pf cx e^^ , this equation determines 



a threshold value e't • If est < Cst', the errors can be 
corrected by successively applying the scheme described 
above. For this kind of concatenated codes, the failure 
probability scales with the number of levels of encoding 
k as 



th 



thN2'= 



(est/e^n 



(E3) 



whereas the number of physical qubits needed to con- 
struct the logical qubits scales as 7*^. In addition to 
the physical qubits needed for construction of the logical 
qubits, a large number of ancillas are typically needed 
to perform the syndrome measurements. These ancillas 
constitute the overhead in the required computational 
resources. 
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Figure 10. Quantum circuit to measure the generators of the Steane code in a traditional architecture that allows only for 
single- and two-qubit gates, and single-qubit measurements. Each of the six generator measurements is realized using four 
CNOT gates with an ancilla, which is in turn encoded using four physical qubits to avoid error propagation. This is the circuit 
we used to compare the error threshold with and without multi-qubit measurements. 



Realization of the Steane code with the ramm 



6, and the circuit contains 6 measurements, we obtain 



In the case of the ramm, the syndromes can be directly 
measured. For simplicity, we assume that one single- 
qubit gate is always performed during the recovery part 
of the circuit. Considering that each qubit is on average 
involved in 24/7 measurements, the total number of time- 
steps required to perform the syndrome measurements is 



i(l) = 6eom , 

,(2) = i.6-5eL = 15e, 



2 



2 

om ' 



-I- 24, 



■ 



(E4a) 
(E4b) 
(E4c) 



These values allow to compute explicitly P(failure, N) for 
the RAMM via Eq. (El). 
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2. Steane's code without multi-qubit measurements 

We want to compare the error threshold in ramm with 
a reference system, where multi-qubit measurements are 
not available. The syndrome measurements are then per- 
formed with the help of ancillas. In particular, the fault- 
tolerant realization of the six syndrome measurements re- 
quires a total of 24 ancillas, each quadruplet being used 
for measuring one of the syndromes [44] (see Fig. 10). 

Each syndrome is measured by first initializing the an- 
cilla quadruplet in a Shor state, which guarantees that 
measuring the four ancillas will not destroy the state en- 
coded in the logical qubit. The second step consists of en- 
coding the syndrome into the quadruplet, which requires 
performing a total of four CNOT gates between different 
ancillas and physical qubits. Since these involve indepen- 
dent qubit pairs, we assume that these four gates are per- 
formed simultaneously. Additionally, we assume that the 
syndrome is measured immediately after the CNOT gates 
and the initialization of the ancilla quadruplet takes place 
already before the syndrome measurements. Because er- 
rors occurring in the ancillas essentially have the same 
effect as the errors in the measurement outcomes, we in- 
clude all possible ancilla errors in Poin(w). 

The initialization of the ancillas to a Shor state is ex- 
plained in Ref. 44. It involves 7 time steps with 5 CNOT 
and 5 Hadamard gates. Moreover, a measurement is re- 
quired to confirm that the Shor state was successfully 
encoded, otherwise the initialization process is repeated. 
We only consider gate and storage errors occurring in the 
initialization of the four ancillas. Each of the ancillas is 
acted on with 13/4 gates on average. 

The syndrome measurements involve 9 time steps and 
each of the 7 physical qubits is acted upon with 38/7 
gates on average, while recovery part only involves one 
single-qubit gate. Finally, we need to take into account 
the errors occurring in any of the 24 ancillas during the 
syndrome block, which contribute to Pom- This way we 
obtain 

Pom (1) = 24 (Pinit + Fsyndromc) , (E5a) 
Pom (2) = 2jLx23 (^p.^.^ Psyndromc)^ i (E5b) 

P,, = f Eg + f e^t + + f est , (E5c) 

with 

^'init = (om + Cdm , (E6a) 

-^syndrome — ~4~^g ^ ~^^st ^ ^om ^ 24 ^g^. (E6b) 

These values allow to compute P(failure, N) in the ab- 
sence of multi-qubit measurements. 

3. Comparison of the error thresholds for the 
quantum memory 

We minimize the probability of failure per time step 
with respect to iV for both implementations of the error 
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1 4 eg/e^t 7 10 

Figure 11. Ratio of the Steane code error thresholds with and 
without multi-qubit measurements as a function of the ratio 
between gate and storage errors, eg /est. The solid, dashed, 
dotted, and dash-dotted curves correspond to ratios eom/est = 
edm/est = 1,2,5, and 10, respectively. 
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Figure 12. Ratio of the computational error thresholds with 
and without multi-qubit measurements as a function of M. 
Here e = eom = edm = eg, with e/tst = 1 (solid), 5 (dashed), 
and 10 (dotted). The range of M starts from 11, because of 
the condition M - No - 1 > 0. 



correction scheme. We characterize the relative proba- 
bilities of errors by fixing the ratios Cg/cst, Cdm/sst and 
Com /fist, and calculate the error threshold for est- Results 
are shown in Fig. 11. We find that for eg — Cdm = ^om = 
est the error threshold of the RAMM is approximately an 
order of magnitude larger than the error threshold of a 
reference architecture that can only perform single- and 
two-qubit operations. The ratio of the error thresholds 
for the different architectures becomes smaller with in- 
creasing measurement errors (larger ratios edm /est and 
eom/est), because it becomes favorable to increase the 
waiting time between the consequent error correction 
steps; but even for eg = edm — eom = lOcst we still find 
that the ramm has an error threshold five times larger 
than the reference architecture. 
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4. Comparison of the error threshold in quantum 
computation 

To estimate the error threshold in quantum computa- 
tion, we consider an algorithm where each qubit partici- 
pates in a two-qubit gate with a randomly chosen other 
qubit after every M > Nq time steps. We assume that 
the syndrome and recovery steps are performed once af- 
ter each two-qubit gate. To estimate the error threshold 
we calculate the probability of failure in any one of the 



which we compute for both architectures using Eqs. (E4), 
(E5). The probability of failure per time step is then 

Pf = P(failure, M) /M, (E8) 

and the threshold for quantum error correction can be 
determined by comparing this probability to the prob- 
ability of failure without error correction. Results are 
shown in Fig. 12. Similarly as in the case of quantum 
memory, we find that the error threshold for performing 
the quantum computation can be an order of magnitude 
larger for the ramm. 

Appendix F: Characteristic energy scales of the 
problem 

We need to satisfy the following chain of inequalities 

Ej^k,hnk,Ag > Ej,o,hno,hwQ > Em > A+, 

Amax > Amin, (Fl) 

where hftk ~ \/^Ej^kEc,k is the plasma frequency of the 
small islands and Ag ~ 100 GHz is the induced gap in 
the nanowire. 

We also have some additional constraints. The charg- 
ing energies and Josephson energies have to be related 
so that Ej/Ec and Ej,o/Ecfl can be varied within the 
range 10 — 50 with the help of magnetic fluxes. By as- 
suming that Ej^/Ecfi = 10 during the measurement, we 
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